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In applying a Ritz modal expansion to the solution of a 
transient response, there is a problem as to how many modes are 
needed to obtain accuracy to within a specified percentage. 

One of us, Chargin, has suggested a method based on the 
characteristics of the forcing function. The method can be 
incorporated into the Ritz generation algorithm such that it will 
automatically monitor, regulate and terminate the process 
according to a specified tolerance. 

FORCING FUNCTION CHARACTERISTICS 

Assume that the forcing function F(x,t) can be repre- 
sented as a product of a spatial function and a temporal func- 
tion; i.e. 

(1) F(x,t ) = F(x) . f (t) . 

Develop a criterion based upon measuring the amount of power 
developed in the forcing function F(x,t). The total power in 
F(x,t) is the product of the "power" in F(x) and f(t) owing to 
the assumption in equation (1). The scheme in outline is to 
measure the temporal power and the spatial power in F(x,t) separ- 
ately then compare the corresponding power in the Ritz modes 
against power in each component of the forcing. Generation of 
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additional Ritz modes will continue until the power criteria are 

met. First a measure is taken of the total power in the temporal 
function. 


Temporal Power 


( 2 ) 


P(x) = 




dt 


where t is the interval over which the transient will act. There 
could very well be a separate temporal function for each spatial 

function. In that case of multiple loadings P(x) of equation (2) 
would be a vector "1" long. 


In order to tailor this power to our use as a guide in 
selecting Ritz modes it will be useful to measure the amount of 
temporal power as a function of frequency. Expand the temporal 

function in a Fourier Series and sum the power versus the expan- 
sion multiple: 


(3) f(t) = a 0 + a ^cos wt + b^sin wt + a 2 cos 2wt +■ b 2 sin 2wt 


cos ut + b sin ut . 
n n 


The power within a band 0 to nu is 
n 

(4) 


P = > fa cos iut + b.sin iwt) = > fa? + b?) . 
n i=0 l 1 1 J j =0 l i ij 


One can compare the amount of power within a given band P with 
the total power P(t). When P^ is within close range of P(x), say 
1%, then the analyst can be satisfied that the frequency range of 
the truncated temporal forcing function is sufficiently broad to 
encompass the temporal requirements of the forcing. 
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P 

For a certain frequency f R = ^ = ^et t ^^- s ^ re ~ 

quency be f . 


Spatial Power 


Now turn to the spatial distribution of the forcing 1 function 
F(x), and develop a measure that is called spatial power. Make 
an additional assumption that all points being loaded have mass. 

(5) Cn(x)3 = CF(x)D T CM3 CF(x) 3 = [ C ia ] [ F ^ } ] ' 


where £c^_J = |f(x)^W| is a 

to be usefual later. H(x) 
loading cases "1" . 


coefficient matrix that will be shown 
is a matrix if there are a number of 
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At this point we have 2 measures of the forcing fuction. 
We have f^ reflecting the desired frequency content and JI(x) 
reflecting the desired level of power to activate the structural 
mass. Now it is time to test the adequacy of the number of Ritz 
modes to repond to the forcing at these power levels. 

First, get a spatial measure of the Ritz modes 0^. A 
simple scheme is to pattern the measure after equation (5) but 
substitute the matrix of k spatial Ritz functions for the post- 
multiply operation instead of the spatial component of loading 
F(x) . 


'6, [' R r\] ■ £ 


CF(x ) 3' 
k=l 


CM3 


N ■t.hJH 


The way to use R^ is to compare each diagonal term of H(x) with 
the corresponding diagonal term of to :ee if the ratio is 
within a specified tolerance; i.e. 


(7) 


IT, 


I 



< e for every 1. 


Keep generating additional Ritz modes 0^ k > r until an r has 
been reached for which every diagonal term satisfies the 
criterion . 


Monitoring Function 

If the 0^ satisfy the spatial power criterion of F(x,t) it does 
not necessarily hold that 0 fc will simultaneously satisfy the 
temporal requirements of F(x,t). Therefore, use the 0^ which 
have met the spatial power requirements and obtain an estimate of 
its frequency content by setting up the frequency equation. Use 
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an estimator that is much less demanding than solving the 
eigenvalue problem. Merely do the following. 

Compute a k— order Ritz generalized mass and stiffness; i.e. 


( 8 ) 




and [kJ = [e^CKI je R 


and construct a test 
threshold frequency, f , 


matrix called 
determined from 


|^S k J which involves the 
the temporal power P . 


(9) [K k - s [sj- 

Decompose and extract the value of the output parameter 

NBRCHG issued by the DECOMP module. Param NBRCHG reports the 

number of negative values on the factor diagonal of which is 

tantamount to the number of sign changes or zero crossings in the 

characteristic equation. If the value of NBRCHG = k, it 

indicates that the frequencies of all k Ritz modes are less than 

f . One would be inclined to want the frequency content of the 
o 

Ritz modes to bracket f ; i.e. some modes with a frequecy ' t Q . 

This implies that one would seek to have the value of NBRCHG to 

be less than the order of matrix S^. This idea can be built into 

the Ritz generation routine as a test as to whether enough modes 

have been generated to within a certain margin ot such that f^ / 

af , where the user specifies a. 
o 


Obviously, one would not want to repeat the eigenvalue 

T 

estimate each time a new Ritz mode is obtained, because e^K 0^ 

and 0^M 0. could become expensive as k becomes large. Therefore 

develop a scheme whereby the spatial power is reduced by some 

factor y > 1 ; i.e. e = e ,./y. Typically, one can use a value 

n o w cd 1 . 0 . 
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of 2 to 4 for y. Once the spatial error e , is satisried, 

H6W 

repeat the eigenvalue estimate test. 

Conclusion 


A scheme hasbeen proposed to monitor the adequacy of a 
set of Ritz modes to represent a solution by comparing the 
quantity generated with certain properties involving the forcing 
function. In so doing an attempt has been made to keep this 
algorithm lean and efficient, so that it will be economical to 
apply. Using this monitoring scheme during Ritz Mode generation 
will automatically ensure that the k Ritz modes 0 that are 
generated are adequate to represent both the spatial and temporal 
behavior of the structure when forced under the given transient 
condition defined by F(x,t). 
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